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We study the Landau-Zener Problem for a decaying two-level-system described by a non-hermitean 
Hamiltonian, depending analytically on time. Use of a super-adiabatic basis allows to calculate the 
non-adiabatic transition probability P in the slow-sweep limit, without specifying the Hamiltonian 
explicitly. It is found that P consists of a "dynamical" and a "geometrical" factors. The former 
is determined by the complex adiabatic eigenvalues E±(t), only, whereas the latter solely requires 
the knowledge of a±(t), the ratio of the components of each of the adiabatic eigenstates. Both 
factors can be split into a universal one, depending only on the complex level crossing points, and 
a nonuniversal one, involving the full time dependence of E±(t). This general result is applied to 
the Akulin-Schleich model where the initial upper level is damped with damping constant 7. For 
analytic power-law sweeps we find that Stiickelberg oscillations of P exist for 7 smaller than a critical 
value 7 C and disappear for 7 > r y c . A physical interpretation of this behavior will be presented by 
use of a damped harmonic oscillator. 

PACS numbers: 34.10.+X, 32.80.Bx, 03.65.Vf 



I. INTRODUCTION 

In many cases one can reduce the quantum behavior 
of a system to that of a two-level system (TLS), which 
corresponds to a (pseudo-)spin one half. The spin-down 
and spin- up state will be denoted by |1) and |2), respec- 
tively. If the TLS is in state |^q) at time to one obtains 
^(t)) by solving the Schrodinger equation 

ihf t \*(t))=H(t)Mt)) (1) 

with initial condition ^(io)) = l^o)- Note, that we allow 
for an explicit time dependence of H. One of the quan- 
tities of particular interest is the survival probability 

P= lim lim |<*(to)|*(i)>| 2 (2) 

t — >oc to — > — 00 

that the system remains in its initial state. For a TLS 
with a level spacing depending linearly on time the result 
for P as function of the sweep rate v has been derived 
approximately by Landau and Stiickelberg Q and rig- 
orously by Zener Q and Majorana Q. P will depend 
sensitively on the t-dependence of H and can not be cal- 
culated analytically, except in limiting cases, only. One of 
them is the adiabatic limit. In that limit it is known that 
\^(t)) converges to a superposition of the adiabatic states 
\ u o,±(t)) which are solutions of the eigenvalue equation: 

H(t)|«o,±(0) = E±(t)\uo,±(t)) (3) 

with E±(t) the adiabatic eigenvalues. Although E + (t) 
and E-(t) may not cross in real time (avoided level- 
crossing) this will happen for complex times t k c , k — 
1,2,. ..,N. 



In case of a real- symmetric Hamiltonian matrix 
(v\H(t)\v') , v,v' — 1,2 which is analytic in t and for 
a single crossing point t c in the upper complex i-plane 
(Irrrt c > 0) it was shown by Dykhne |5| (see also earlier 
work by Pokrovskii et al. jgj) that 

P = cxp[-2 Im z(t c )} (4) 

in the adiabatic limit. The new variable z(t) is given by 
t 

z(t) = J dt'[E + (t') - E_(t% (5) 
'0 

Davis and Pechukas 7] have performed an exact proof 
of result (0}, JSJ. Particularly, these authors have proven 
that the pre-exponential factor equals one. Therefore it 
is sometimes called the Dykhne-Davis-Pechukas (DDP) 
formula. For more than one crossing point with lmz k = 
ImzCtc) > a generalization of (0} has been suggested 
0, and tested by Suominen and co-workers (Ref. 
and references where- in) . A rigorous prove of the general- 
ization of DDP-formula including even hermitean Hamil- 
tonians has been provided by Joye et al. . More than 
one crossing point leads to interferences which generate 
oscillations in P as function of control parameters, like 
the sweeping rate (see below). 

For Hamiltonian matrices which are not real- 
symmetric, but hermitean, Berry ^lj and Joye et al. 
made an interesting observation which is that P obtains 
also a "geometrical" factor besides the "dynamical" one, 
Eq. QJ, where the former also depends on the crossing 
points t k c , only. For those who are less familiar with this 
kind of physics let us explain the choice of this nomen- 
clature. Below we will see that one of the factors of P 
is entirely determined by the adiabatic eigenvalues and 
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the other by the adiabatic eigenstates. Since the former 
is important for the time evolution it is called "dynami- 
cal" whereas the latter is related to the geometry in the 
Hilbert space, particularly through a condition for par- 
allel transport fEa. I21f> . and accordingly it is called "ge- 
ometrical" . 

TLS will be influenced by their environment, e.g. by 
phonons. The spin-phonon coupling leads to dissipa- 
tion of the (pseudo-)spin dynamics which will influence 
the probability P. Although there exist microscopic 
models for the spin-boson system |13(, and simplified 
models where the bath is described by fluctuating fields 
[lH HH llH Il7l Il8| , we will use a dissipative Schrodinger 
equation. This will be achieved by using a non-hermitean 
Hamiltonian for the TLS. A particular version of such a 
model has been suggested by Akulin and Schleich |l9j . In 
their model, called AS-model in the following, the upper 
level (at the initial time to) experiences a damping (see 
section III). 

The survival and transition probability for non- 
hermitean TLS-Hamiltonians has already been investi- 
gated by Moyer j2jj. This has been done by mapping 
the original differential equation to the Weber equation, 
which can be solved exactly. By use of the Weber equa- 
tion as the appropriate "comparison equation" it was 
shown how the DDP-formula, Eqs. Q and 0, can be 
extended [2(j. However, this extension does not con- 
tain a "geometric" contribution, although one expects 
that it exists similarly to what has been proven for her- 
mitean matrices |ll|, [l2|- On the other hand Garrison 
and Wright [2l| have investigated the geometrical phase 
for dissipative systems but not the non-adiabatic transi- 
tion probability. 

It is one of our main goals to derive a generalized 
DDP-formula in the adiabatic limit containing a "geo- 
metrical" and a "dynamical" contribution for a general 
non-hermitean TLS-Hamiltonian. We will demonstrate 
that both contributions consist of a universal and a non- 
universal part. The former depends only on the complex 
crossing points whereas the latter requires the knowledge 
of the complete time dependence of H . Instead of us- 
ing a "comparison equation" we apply the concept of a 
superadiabatic basis, put forward by Berry |22|. to non- 
hermitean TLS-Hamiltonians. As a result we will find 
that the "dynamical" contribution to the non-adiabatic 
transition probability (which equals the survival proba- 
bility in the adiabatic limit) is determined by the com- 
plex, adiabatic eigenvalues E± (t) , only. The correspond- 
ing "geometrical" part solely requires the knowledge of 
a ±(0)> the ratio of the components of each of the adia- 
batic eigenstates. 

A second motivation is the application of our results to 
the AS-model. It has been shown that the survival prob- 
ability P does not depend on the damping coefficient 7 
of the upper level, provided the bias of the TLS varies 
linearly in time, and the coupling A between both lev- 
els is time-independent [l9| . Therefore it is interesting 
to investigate non-linear time dependence and to check 



whether or not P remains insensitive on 7. For non- 
linear time dependence more than one complex crossing 
points may occur, such that interference effects can gov- 
ern the dependence of P on the sweeping rate jlQj. Spe- 
cific examples with 7 = for which this happens were 
discussed in recent years [jj |2^| . There it was found that 
critical values for the sweeping rate exists at which the 
survival probability vanishes, i.e. complete transitions 
occur between both quantum levels. Consequently, one 
may ask: Are these complete transitions reduced or even 
suppressed in the presence of damping? 

Our paper is organized as follows. The next section 
will contain the general treatment of the non-hermitean 
Hamiltonian and the presentation of the generalized 
DDP-formula. In section III we will apply the results 
from the second section to the AS-model with power 
law time dependence. The results for the AS-model for 
power law sweeps can be interpreted by the dynamics 
of a damped harmonic oscillator. This will be shown in 
section IV. A short summary and some conclusions are 
given in the final section. 

II. GENERAL FORMULA FOR 
NONADIABATIC TRANSITION PROBABILITY 

In this section, we will derive a generalized DDP- 
formula for the non-adiabatic transition problem of a 
decaying TLS. The Hamiltonian can be represented as 
follows 



H{St/K) = -Y,B j {5t/%)o- j (6) 
3=1 

with <jj, the Pauli-matrices and Bj a time dependent 
field. 5 > is the adiabaticity parameter. Because this 
model should be dissipative, at least one of the Bj must 
contain a nonzero imaginary part. Accordingly H is non- 
hermitean. In the following we will assume that Bj is 
analytic in t. Introducing a new time variable 

r = 5t/h (7) 

Eq. Q becomes 

i5d T MT)) = H(T)\*(T)), (8) 

where d T = d/dr. |^(r)) can be expanded with respect 
to \u) 

2 

iw>-E c ^ r )H- ( 9 ) 

With |\&(to)), the initial state, its survival probability is 
P = P{5)= Urn lim |(*(r )|*(r))| 2 . (10) 
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Note that P is the survival probability with respect to 
the diabatic basis. With respect to the adiabatic basis P 
is the nonadiabatic transition probability. 

To calculate P for 5 -C 1 we introduce the adiabatic 
basis of H(t). This can be done as in Ref. |2l| where a 
biorthonormal set of right-eigenstates was used or alter- 
natively by use of left- and right-eigenstates. We will use 
the latter, as it turns out to be more elegant. Let 



(ii) 



be the adiabatic rig/if-eigenstates. They are solutions of 



H(t)\u ,±(t)) = E ± (t)\u ,±(t)) 



(12) 



with E±(t), the adiabatic eigenvalues. Note that E±(t) 
are complex in general and that the norm of \uq,±(t)) and 
of ^(t)) is not conserved, since H(r) is nonhermitean. 
Following Berry |22| , we introduce a hierarchy of supera- 
diabatic right-eigenstates \u nj ±(r)) , n — 0,1,2,... and 
expand the solutions \^±(t)) of Eq. JSJ with respect to 
the superadiabatic basis: 



|*±(t)> =exp 



dT'E±{T r ) 



m—0 



Substituting \ 1 &±(t)) into Eq. JHJ yields the recursion 
relations 



[H(T)-E a (r)}\uoAr))=0 



a = ± 



(14) 



%d T \u m ^^(r)) = [H(t) - E a {T)}\u m . a {T)) , m> 1. 

(15) 

Eq. {H} is already fulfilled, due to Eq. {TJJl. To make 
progress we introduce the adiabatic Ze/i(-eigenstates 



(16) 



which are solutions of 

(Zo,„(t)\H(t) = E a (T)(u Q , a (T)\ (17) 
and are normalized such that: 

(wo,o-(r)|Mo, CT '(r)) = 5 aa <. (18) 

Let be: 



, v 4(r) 



(19) 



the ratio of the components of the adiabatic right- 
eigenstate \uq^(t)). Then it is straightforward to prove 
thatmm 



-ct T (r), v = 1 
[a + (r)-a_(r)]ei(r) \ L ^ = 2 



(20) 



which defines the left-eigenstate from the right- 
eigenstate. Multiplication of Eq. I|15|) for m = 1 with 
(wq, ct (t)| leads to 



(uo,cr(T)\d T \uo,cr(T)) = , (7 = ±. 



(21) 



This is the condition for "parallel transport" jl lll2 1| now 
generalized to nonhermitean Hamiltonians. 

In order to solve recursion l|15(l we expand \u mja {r)) , 
m > 1 with respect to \uq i<7 (t)): 

\umAr)) = a* (r)|uo,-(r)) + ^(r)K+(r)). (22) 

Substitution of Eq. l(2"2"|) into Eq. lfT5|l and multiplying by 
(wo, ct (t)| yields with Eqs. (HU), JTSJ for m > 1: 



a m -i( T ) = -K-(r)6 m _i(r) 



(23) 



&m-l( T ) = - K +( T K-lW - - E_( T )]b m (T), 

(24) 

where ' denotes derivative with respect to r and 



r (r) = (t2o, (T (T)|(i T |uo,_ (T (T)) 



(25) 



are the nonadiabatic coupling functions, responsible for 
the nonadiabatic transitions. If k<j(t) = 0, we get from 
Eqs. (123 and (El 



a m (T~) = a rn (T ) , 6 m (r) 



S+(r)-£_(r) 



&m-l( T )- 



(26) 

Similar equations follow for a+ (r), &^(r), which however, 
will not be needed. Next we fix the initial condition for 
IMt)>: 



|# ct (t )) = |n 0i<7 (T )), 



(27) 



i.e., we start in the adiabatic right-eigenstates. From 
Eqs. (JT3J, (1221 we find immediately for a = — 



a (r) 
a m( r o) = 



1, Mt)^0 

6-(t ) = 0, m>l. (28) 



such that Eq. (|2l)|l implies a rn (r) = 0, 6~ (r) = 0, m > 1 
provided «±(t) = 0. This makes obvious the absence of 
nonadiabatic transitions. 

The next step is the calculation of k ct (t). For this we 
need e^(r), which can be determined from J2J. As a 
result we find 



4(t) = 4(t ) 



exp 



dr' 



q.'+(t') — a_(r') 



, (29) 



and taking Eq. (|20fl into account we obtain the general 
result 



K<t(t) 



a—r- — -exp 

4( T o) 

a + (r) - a_(r) 



dr 



, A + (^)+a-(0 
q; + (t') — a_(r') 



exp 



a + {T') — a_(r') 
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where the expression has been split into a r-independent 
(first line) and a r-dependent factor (second line). Fol- 
lowing Berry we truncate the series, Eq. JT^Jl, at the 
n-th level 



l*n, CT (V)) = exp 



/ dr'E^r') 



(31) 



m=0 



and expand |\I>(t)): 

|*(t)> = E c n Ar)\*nAT))- (32) 

<T = ± 

As initial condition we choose: 

|*(t )) = |*_(r )), (33) 
which is equivalent to 

c n _(r ) = 1, <V + (r ) = 0, n -> oo. (34) 
Introducing a corresponding truncated state 

n 

(*„,*(t)| - /»,„(r) E <H<wM|, (35) 

m=0 

where the r-dependent prefactor f n ,a{A has not to be 
specified we obtain an equation of motion for c Ili(7 (r), 
after Eq. (|52)l has been substituted into Eq. ©: 



i<5c„ i<T (r) = E H n . aa > (r)c„ >CT / (r) 



(36) 



with 



rx"=± 

n n .^(T) = (y nia (T)\H(T)-i5d T \y n , a ,(T)}. (37) 
Eq. (|36|) can be rewritten as an integral equation 



Cn,a(r) = c n)O -(r ) + ^ / dr ' ' H rL . <J - <J {T l )c n ^ (T { T ') 



x exp 



H n . aa (T ) 



(38) 



Apart from the truncation, Eq. I|31|) . the results are still 
exact. Eq. (|38|l simplifies in the adiabatic limit 6 — ► 0. 
In leading order in 5 we get from Eqs. (|18|l . I|31() and 1|35|) 



^„ ;o - CT '(r) = /n, CT (r) exp 



~ / dr'S ff (r') 



(39) 



[if(r) - i<ta T ]|*„, CT /(r)) can be found in Ref. HJ. Mul- 
tiplying by (^ n)<T (r)| and making use of Eqs. (fT2jl . iJTSJl, 
(J22J and leads to 

Un,aa\T) = -F* 1 [E a {f) " ^ (/)] 



x exp 



T 

-j | dr%,(r') 



x < +1 (r)<5 ff ,_+J£ +1 (r)5 CT)+ +G(r^) (40) 
from which follows 



ff„;^(r)--<5" +1 
x [£ CT (r)-£ CT ,(r)] 



< +1 (t)<5 ct _ + 6^+1^, 



exp i 1 |dr' [^(r') - £ a ,(r')] 1 + 0(<T +2 )(.41) 



Note that the prefactor f n ,A T ) has cancelled. The diag- 
onal elements of H„(r) are of order S n+2 and the non- 
diagonal ones of order S n+1 . Therefore it follows from 
Eqs. (ESJ), (HP, (Ell) and JSHJl 



/ \ ^ .,„el_(ro) 
Cn,+ (Tj = *<3 -f7 — rexp 



dr 



, a + (r')+A_(^) 
a+(r') — a_(r') 



X / ^d- +1 (,0 [a+(r 2 + ^- (T/)] [^(rO-^(rO]. 
■ a+(r ) + a_(r ) 



x exp 



a+(r") - a-(r") 



xexp<( - / dr [£+(r )-£7_(t )] 



(42) 



The time dependence of H(t) is chosen such that 

lim |w ,-(to)) = |1) , lim uo,+ (t )) = |2) 

TO — * — OO TO — * — OO 

lim |«o,-(t)> ~ |2) , lim |n 0)+ (r)) ~ (43) 

r — >oo r — >oo 

Note that the adiabatic states at initial time To are nor- 
malized. Since Eqs. J23, lO and g3J imply 



lim |*( T0 )) = |1) 



(44) 



we obtain from Eq. (|10|l for the nonadiabatic transition 
probability in leading order in S 



c n , + (oo)(l|u ,+ (oo)) exp 



dT'E + ( T r ) 



(45) 
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where we used (l|ito._(oo)) = 0, due to Eq. I|43|l . Substi- 
tuting c„ !+ (oo) from Eq. with t = — oo into Eq. l|43|) 
we get with Eqs. (tTT|l . and linv^-oo el (r ) = 1 
(due to Eq. 103} ) 



P(5) = cxp 



pus , * pn. 



I dra- +1 {r) a+{T) ^ [E + (r) - E_(r)} 
a+(T) 

,a+(r')+a-(r') 



x exp 



dr 



aq_(r') — a_(r') 



xexp<(- I dr'[E + {r')-E^T') 



(46) 



with the nonsingular "geometrical" and "dynamical" 
contribution 



F" s = 2Rc 



d+(r) 



a+(r) — a_ (r) 



d_(r) 



a+(r) - a_(r) 



(47) 



and 



-21m 



OO U 

y drE+ir) + J drE^{r) 



(48) 



respectively. Note that FJ 18 = 0, for a hermitean Hamil- 
tonian, since E a (r) are real. The expressions for F™ s and 
F£ s put some constraints on if (t), because both quanti- 
ties should be larger or equal to a constant c > — oo, 
which requires that \thE±{t) decays fast enough for 
r — > ±oo. 

The r-integral in Eq. I|4tj[l is dominated by the sin- 
gularities of E + (t) — £L(r), for 6 — > 0. The adiabatic 
eigenvalues and a±(r) have the form 



1 



B±(r) = - T(r)± Vr 2 (r) - 4D(r) 



(49) 



/ ^ -gn(r) + if 22 (r) ± y/r 2 (r) - 4D(r) 

where T and D, is respectively, the trace and the deter- 
minant of the Hamiltonian matrix H vu i = (y\H\v'). Ac- 
cordingly, the singularities are the branch points r c (fc), 
k = 1, 2 , ... of E + (t) — E_(t). Introducing a new vari- 
able fill El 



r 

r(r)= | dr'[F + (r')-F_(r')] 



(51) 




us 



/ 



FIG. 1: it-dependence of the adiabatic eigenvalues for a power 
law sweep w(u) = u 3 and 7 = 0.5 < ■y c — 1 (a) real part of 
E±(u), and (b) imaginary part of E±(u) 



it is shown in the Appendix A that after taking the limit 
n — > 00 the nonadiabatic transition probability is given 
by 

2 



P(S) = exp 



ns 1 77ms \ 



expire**' 



(fc) 



with the singular "geometrical" contribution 



2 c (fe) 

F°(k)= f dz 



dz 



■{*) 



da.— 
dz 



14.(2) - a-(z) 



(52) 



(53) 



and the singular points z c (k) — z(r c (k)), which are above 
the contour C — {z(t) \ — 00 < r < 00}. The final result 
of this section, Eq. I|52(l . is the generalization of the DDP 
formula (as it has been rigorously proven for hermitean 
TLS-Hamiltonians [ToT ] ) to nonhermitean ones, describing 
dissipative TLS. The reader should note that the use of 
the superadiabatic basis leads to a pre-exponential factor 
in Eq. (|52|l which is equal to one, which is identical to the 
case without dissipation. The result, Eq. I|52|). exhibits 
that the "dynamical" contributions follow from the adi- 
abatic eigenvalues and their branch points, whereas the 
"geometrical" contributions involve u±(t), only. If we 
parametrize for a TLS with hermitean Hamiltonian the 
external field components Bj, Eq. Jf)J), as it has been done 

and that Eq. JSJ 



in Ref. [TJ 

becomes: 



one recovers that F™ s 



r c {k) 

F° = J dr 0(t)cos9(t) 


in agreement with the result in Ref. |ll| . 

III. APPLICATION TO THE 
AKULIN-SCHLEICH MODEL 

The AS-model is given by [H| 
1 



(54) 



H(t) 



[W{t)a z + Aa x + i~f(a z - c-a) 



(55) 
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FIG. 2: Same as Fig. 1, but for 7 = 1.1 > 7 C 



with the external field W(t), the tunnelling matrix el- 
ement A, and the damping constant 7 > of level 
|2)=| f)< Co is the 2x2 unit matrix. Let us introduce 
dimensionless quantities: 



w(u) 



W(t) 



vt 
A' 



7 



2 

A' 



A 2 

Tiv 



(56) 



Note that the time variable r of the previous section 
is not dimensionless. After the replacement of t by it, 
Eq. £Q) takes the form of Eq. © with: 



From Eqs. (|49Jl and l|5()(l it follows immediately 



E±{u) = - —vy ± y (w(u) + ry) 2 + 1 



(57) 



(58) 



a±(u) = -(w(u) + ij) ± ^/{w{u) + n) 2 + 1, (59) 

where the branch of the square root has been chosen such 
that > 0, for x > 0. Fig. CJi and Fig. exhibit 
Rei?±(u), and ImE±(u), respectively, for an analytical 
power law sweep w(u) = it 3 and for 7 < 1. The corre- 
sponding result for 7 > 1 is shown in Figures^ and [3). 
In the following we will consider crossing sweeps, only. 
For those it is 



lim w(u) = ±00. 

u — ►ioo 



(60) 




FIG. 3: Branch points (open circles) M^(fc) of E + (u) — E^(u) 



for a power law sweep w(u) 



and 7 < 7 C . The thick solid 



lines are the branch cuts. The radius of the inner and outer 
circle is (1 — 7) 1 / 3 and (1 + y) 1 ^ 3 , respectively 



is fulfilled if w(u) decays faster than u -1 / 3 . Otherwise 
F2 S — 00 which makes P(8) to vanish. Figs. ^ an d [21 
demonstrate that there exists a critical value for j c . For 
< 7 < 7 C — 1 we have Re[E + (u) - E_(u)] > for all u 
and ImEuiu) is continuous whereas Re[£"+(ii) — E—(u)] 
vanishes if w(u) = and lmE a (u) becomes discontinu- 
ous, provided E±(u) are defined by Eq. (|58l . 

The nonsingular geometrical part, Eq. (|47|l . can be cal- 
culated without specifying the u-dependence of w. Sub- 
stituting a±(u) and a±(u) from Eq. (|59|l into Eq. I|47|l . 
both integrals in Eq. 147|l . become a sum of two integrals. 
One of them can be calculated by the introduction of a 
new integration variable £ = w + 27 and the other by 
noticing that its integrand can be rewritten as a deriva- 
tive of a logarithm with respect to u. Without restrict- 
ing generality we assume that w(0) — 0. Then we obtain 
with Eq. pi 



Returning sweeps for which lmi u ^± 00 w(u) = —00 (or 
+00) can be treated analogously. It is easy to prove that 



E±(u) 



w(u) + 



1 



If 



2w(u) 2w 2 (u) 



0(vj- 6 {u)) 



for u — > ±00 and 



E±(u) 



-w(u) — 2ij 



2w(u) 



0{w-\u)) 



(61) 



(62) 



for u — > =Foo. 

From Eqs. (JHU and igHjl it follows that F 7 d ls is fi- 



nite provided J du'w 3 (it') exists for 



±00. This 



f; s ( 7 ) = 2Re 



in(»7 + vi-r) 



(63) 



The nonsingular "dynamical" and both singular con- 
tributions require the explicit it-dependence of w. As 
said above we will consider crossing sweeps only. There- 
fore we restrict ourselves to power law sweeps w{u) = u n 
with n > and n odd. n should not be confused with 
the truncation number n in the previous section. Since 
w(— u) = —w(u) we can rewrite F£ 3 as follows: 



F d " s ( 7 ) = 2 du 7 - I mv /(w(u) + i7) 2 + l . (64) 



It is easy to see that 

F™(0) = 0, *T(G) = 0, 



(65) 



for 7 = 0. Hence, the nonsingular contributions to the 
nonadiabatic transition probability vanish if there is no 
dissipation. In this case the result 1|52|) reduces to that 
found by Berry for hermitean Hamiltonians and for 
a single complex crossing point contributing to Eq. (|52[1 . 
What remains is the determination of the singular points 
u c (k), k = 1, 2, . . . and the calculation of z c (k) and Fg{k). 
These singular points are the branch points of E + (u) — 
E- (u) . Their location depends on whether < 7 < 7 C or 
7 > 7 C = 1. Let us start with the first case < 7 < j c . 
From (u n + ij) 2 + 1 = 0, n odd, we find 




1 5 10 15 20 25 




5 10 15 20 25 



: (fc) =±(l=F7) 1/,r exp 





, 27rV 




Kfc— 

n J 



(66) 



for k = 0, 1, . . . , n — 1, which are shown together with the 
branch cuts in Figure 3 for n = 3. From Eqs. I|51|) and 
(|66|l we obtain the corresponding singular points in the 
complex z-plane: 



z t 0) = ±^n(7)exp 



7T , 2TT 

h k — 

2n n 



where 



(1=F7) 



si 1 /* 



»±(7) = J dxy/l-tf±x")*. 



(67) 



(68) 



Since the mapping z(u) is analytic in the complex u- 
plane, except at the branch lines, it is conformal. Accord- 
ingly, for those uf(k) which are in the upper w-plane the 
corresponding z^(fc) will be above the integration con- 
tour C and therefore will contribute to P (see end of the 
second section). After the determination of the singular 
points we can proceed to calculate their "geometrical" 
and "dynamical" contribution to P. From Eqs. i|53|) and 
it follows: 



u±(fe) 



W) = 



(III 



a+(u) — a.-(u) 



a+(u) — a—(u) 

w(uf (fc))+i7 

dC- 



1 



1 



■ ln(-a_ (u^ (A:))) +ln(i7^ 



Vl- 7 2 )(69) 



Because (w(tt^(fc)) + 17) + 1 = we get from Eq. 
that a-(u^(fc)) = 1 such that 



Fail)- 



(70) 



F°(k) =ln(*7+ yi-7^)+*7ri 

The reader should note that i 7 ^ is independent on fc. 
Consequently it can be taken in front of the sum in 



FIG. 4: Comparison of the numerical exact (solid line) and 
the asymptotic result, Eq. I|71^. (dashed line) for P(e, 7) and 
a power law sweep w(u) — u 3 . (a) 7 = 0.1, (b) 7 = 0.3, (c) 
7 = 0.7 and (d) 7 = 0.9 



Eq. 152|) which yields exp(2ReF* ) and just cancels the 
non-singular "geometrical" factor exp(— 2ReF™ s ), due to 
Eq. I|63|) . Therefore we find that no "geometrical" factor 
occurs for the AS-model. This will change if we apply an 
additional time- dependent field in the x- and y-direction. 
What remains is the calculation of the singular "dynam- 
ical" factor. Because we are interested in the adiabatic 
limit 5 — > 0, we have to take into account in Eq. I|52() 
those singularities in the upper z-plane with smallest 
imaginary part. These are zj (k = 0) and zj(k = 
(n - l)/2), for which Rez+(0) = -Rcz+((n - l)/2) and 
(of course) Imz+(0) = lmz+((n - l)/2). Using Eq. l(57|l 
with fc = and k = (n — l)/2 we obtain finally: 



P 



4 cos 2 (eft+(7) 
x exp [— 2e/i+(7)l sin 



cos— lexp[-eFr(7)] 

7T 

2n' 



(71) 



Let us consider linear sweeps, i.e. n = 1. Then there 
exists only one singularity it+(0) =1(1 — 7) in the upper 
-u-plane and Eq. lj5"2Tl reduces to 



exp 



-i[Fr + 2Imz+(0) 



(72) 



The exponent can be calculated analytically by using 
u + ¥y as an integration variable. As a consequence one 
finds that the 7-dependence drops out from the exponent. 
With 5 = eT 1 one obtains 



P^ e" 



7T_ 7rA 2 

2 C _ 2?iu ' 



(73) 



consistent with the finding in Ref . [ljj . In order to check 
the validity of Eq. O, we have solved numerically the 
time dependent Schrodinger equation in order to deter- 
mine P. A comparison between the numerically exact 
and the asymptotic result, Eq. (|71|l . is shown in Fig. 0] 
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u;(D 




FIG. 5: Same as Fig. 3 but for 7 > 7 C and without branch 
cuts. The radius of the inner and outer circle is (7 — l) 1 / 3 
and (7+ 1) 1/3 . 




15 10 IS 20 25 




15 10 15 20 25 




FIG. 6: Same as Fig. 4, however for 7 > 7 C and with the 
asymptotic result, Eq. 1751 . (a) 7 = 1.01, (b)7 = 1.1, (c) 
7 = 50 and (d) 7 = 500 



for n — 3 and four different 7- values. We observe that 
the deviation between both results, e.g., for e = 5 and 
7 = 0.3, is about 1.6 per cent, only. Similarly good 
agreement has been found for n > 3. From Eq. I|71|) 
it follows that there exist an infinite number of critical 
values (7)1 v = 1)2,3,... at which the oscillatory 
prefactor in Eq. I|71l) vanishes. From this we can con- 
clude that these Stiickelberg oscillations proven to exist 
for TLS without dissipation [Tjj an d discussed later in 
Refs. 0, for 7 = survive even in presence of dissi- 
pation, provided 7 < 7 C = 1. Indeed, we will see below 
that they disappear for 7 > j c . It is not only the survival 



i(T 
10"' 



10" 



10 



10 



critical damping 1 _^ 






* 

\ / — "S P s 

*v ^ V V 1/ / 

- ~- — \/<\/ r w/ / 

~^~>>\f * VI/// 




- 


■ \/' "\ V \f\n i / / 
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1 1 

i 







0.1 



FIG. 7: 7-dependence of P(e, 7) for w(u) — it 3 and e = 
10,15,20,25,30 (from top to bottom). Numerical exact re- 
sult (solid line) and the asymptotic one (dashed line) 



of the oscillations, but also the survival of the complete 
transitions from state |1) = | 4) to state |2) = | f) found in 

Refs. mm for 7 = 0, as long as 7 < 7 C . 

Now, we turn to the second case 7 > 7 C . For this case 
we find: 



: (fc) = -(7=F 1) "exp 





, 27rV 









(74) 



for k = 0, 1, ... ,n — 1, which are shown in Fig. [S] for 
n = 3. The main difference to the case < 7 < 7 C 
is that there is exactly one singular point among uf(k) 
denoted by vP c for which zj? = z{vP c ) is on the real axis in 
the complex z-plane. Using the definition, Eq. I|58l) . of 
E±(u), Fig. |21 demonstrates that E±(u) is discontinuous 
on the real u axis. There seem to exist two possibilities 
to deal with this problem. First, after having chosen the 
branch cuts in the complex it-plane one has to deform 
the integration contour along the real w-axis sucht that 
u = is above that contour and that no branch cut is 
crossed. This kind of reasoning was used by Moyer in 
Ref. [2(j. Second, one could define E±(u) such that they 
are analytic in a strip around the real u-axis. This can 
be done by interchanging E + (u) and E^(u) for u < 0. 
This has the consequence that the contour z(u) for u real 
is in the right complex z-plane, starting e.g. above the 
positive real axis for u = — 00 , going through z = for 
u = and then continuing below the positive real axis 
for u — > +00. This contour would enclose zq if Rezo > 0. 
Whether it can be closed such that the closure does not 
make a contribution is not obvious. Since we are not 
sure how to solve this problem in a rigorous manner, we 
have assumed that zq is the leading contribution to P, 
Eq. for S 0. Since \e lz °/ 5 \ = 1 and due to the 

absence of a "geometrical" contribution we obtain: 



P = exp[-6Fr(7)], 



(75) 



with F^ s (7) given by Eq. (|64|l . A comparison be- 
tween the e-dependence of the numerically exact and the 
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(a) 




(b) 




Comparison of the critical values (7) f° r ( a 



The numerically exact re- 



FIG. 

w(u) — u b and (b) w(u) 
suit is shown by the open circles and the result obtained for 
the oscillator model is depicted by the full circles. The solid 
lines are a guide for the eye. 



asymptotic result, Eq. J7SJ), is presented in Fig. Again 
we find a very good agreement already for e > 1. This 
strongly supports the correctness of our assumption that 
zo is the most important singularity. Eq. I|75|l reveals that 
the Stiickelberg oscillations as function of e have disap- 
peared. We stress that both asymptotic results, Eq. l|7T|) 
and H75|) , are valid for all 7 with < 7 < 7 C and for all 7 
larger than j c , respectively, provided e is large enough. 
This is demonstrated in Fig. 0for different e. 



IV. INTERPRETATION BY A DAMPED 
HARMONIC OSCILLATOR 

In this section we will give an intuitive explanation 
of the Stiickelberg oscillations and will present an ap- 
proximate calculation for the critical values £^(7) for 
power law crossing sweeps w(u) = u", n odd. Close to 



the resonance at u = we may neglect w(u). Then the 
time dependent Schrodinger equation for the amplitude 
of state II) 



i / du'w(u') 



becomes 



with: 



ci(u) = ci(u) exp 



Ci + 2fj,cx(u) + lo ci(u) = 



A* = 



n 



(76) 



(77) 



(78) 



Let ttrans be the Landau-Zcncr transition time. In 
the adiabatic limit it is well-known that ttrans = A/u. 
Eq. H56J) yields ittrans = 1< Therefore we will require as 
initial conditions: 



Cl (-litre 
Cl(— litre 



(79) 



Eq. (|77|) is the equation of motion for a damped har- 
monic oscillator which can easily be solved. The special 
solutions are exp [iui± (e, 7)11] with 



w(e,7) 



1-7 2 



(80) 



This result makes obvious the existence of a critical 
damping j c = 1. For < 7 < j c and 7 < 7 C the oscilla- 
tor is underdamped and overdamped, respectively. This 
qualitative different behavior is the origin of the different 
e-dependence of P for < 7 < 7 C and 7 < j c , found 
in the third section. This relationship can be deepened 
more by calculating £^(7). Having solved Eq. (|77l) with 
initial conditions, Eq. I|79(l we approximate P by: 



P = |ci(+u trs 



-II 



|ci(+lt tr 



+1)| 2 . (81) 



The zeros (with respect to e) of P yield e.£' (7). A numeri- 
cal solution of the corresponding transcendental equation 
leads to the results shown in Figure 8 for w(u) = u n with 
n = 5 and n = 51 and 7 < j c . Figure 8 also contains the 
result from a numerically exact solution of the time de- 
pendent Schrodinger equation. Comparing both results 
we observe that the agreement for n = 5 is qualitatively 
good, but quantitatively less satisfactory. However, in- 
creasing n more and more leads even to a rather good 
quantitative agreement, as can be seen for n — 51. This 
behavior is easily understood, since w(u) within the tran- 
sition range (—1,1) becomes practically zero for n large 
enough. Figure 8 also demonstrates that increases 
monotonically with 7 which is related to the decrease of 
Reu;(e,7) for increasing 7. The oscillator model can also 
be used to determine a lower bound for e < ^\'j — 0). For 
"trans = 1 one gets 



(7 = 0) > 



(82) 
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such that 6^(7) > c { c\l) > c { cHj = 0) > tt/2, for all 
7. It is interesting that the lower bound for e is 
similar to that obtained from the inverse Landau-Zencr 
problem |25| . There, the t-dependent survival probabil- 
ity P(t; e) is given and W(t;e) is determined analytically 
from P(t;e). If P(t;e) = P(t;eu), with u and e from 
Eq. I|56|) . varies from one (for t = —00) to zero (for 
t = +00), it is found that a solution W(t; e) of the inverse 
problem only exists, if 

e > 1. (83) 

The latter inequality, as well as inequality (|82H implies 
that the ratio itransAtunnei of the transit time itrans = 
A/v and the time period of coherent tunneling itunnei — 
K/v, which equals e, is of order one. It is obvious that 
complete transitions can not occur if ttrans is too small 
compared to itunnei, i.e. for e <C 1. In that case the 
quantum system does not have time enough to tunnel 
from the initial state |1) to state |2). 

V. SUMMARY AND CONCLUSIONS 

Our main focus has been on the derivation of the non- 
adiabatic transition probability P(e) for a dissipative 
two-level system modelled by a general non-hermitean 
Hamiltonian, depending analytically on time. Follow- 
ing for the hermitean case Berry's approach by use of a 
superadiabatic basis we have found a generalization of 
the DDP-formula. Besides a "geometrical" and a "dy- 
namical" factor, completely determined by the crossing 
points in the complex time plane, we also have found 
a non-universal "geometrical" and "dynamical" contri- 
bution to P. The latter require the knowledge of the 
Hamiltonian's full time dependence and are identical to 
zero in the absence of dissipation. Without specification 
of the TLS-Hamiltonian, we have shown that both "ge- 
ometrical" contributions can be expressed by a±(u), the 
ratio of the components of each adiabatic states \uq 7 ±(u)) 
in the basis v = 1, 2, and both "dynamical" ones by 
the adiabatic eigenvalues E± (u) , only. In this respect our 
result for P(e) is independent of a special parametriza- 
tion of the Hamiltonian matrix. Although the result in 
Ref. is not in such an explicit form like Eq. (|46|l the 
existence of this nonsingular "dynamical" contribution 
has already been stated there. However, the nonsingular 
"geometrical" part, Eq. I|47|l . has not been found in that 
paper. 

As a physical application we have studied the AS-model 
[T^| . This model describes a dissipative TLS where the 
initial upper level is damped. In 19] it has been shown 
that the probability P for a linear time dependence of the 
bias does not depend on the damping constant 7 for all e. 
Our results demonstrate that this is not generic. For in- 
stance, nonlinear power law crossing sweeps generate a 7- 
dependence of P. For such sweeps a critical value j c = 1 
exists. Below 7 C the non-adiabatic transition probability 
oscillates and vanishes at critical values fig (7), and for 



7 > 7 C the oscillations are absent. Hence, the existence 
of complete transitions at an infinite set of critical sweep 
rates still holds for all 7 below j c . In the section IV we 
have shown how the oscillations and their disappearance 
for 7 > 7 C can be qualitatively explained by a damped 
harmonic oscillator. For power law sweeps with rather 
large exponent, e.g. n — 50, this description becomes 
even quantitatively correct. No doubt, it would be inter- 
esting to study a microscopic model of a TLS coupled to 
phonons, e.g. a spin-boson-Hamiltonian as in Ref. [Isl ]. 
in order to check whether the e-dependence of P exhibits 
oscillations for power law sweep with n > 1 and small 
enough spin-phonon coupling. Another question con- 
cerns the interaction between the TLS which have been 
completely neglected in our present work. That they can 
play a crucial role was shown recently |26( . Whether the 
oscillations still exist in the presence of interactions be- 
tween the TLS is not obvious. 
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Appendix 

In this appendix we will describe how the asymp- 
totic result, Eq. has been derived from 1|46[1 . Al- 
though we follow Berry's approach 0, [2^ we repeat the 
most important steps since the non-hermitean property 
of H does not allow the simple parametrization used in 
Ref. [13 and is not of the form of Eq. @ or Eq. (gZJ 
of Ref. [23. Nevertheless we will recover the same 
universal recursion relation for the coefficients a~ L (r) as 
found by Berry J^]. In order to show how Eq. I|52|) 
can be obtained from Eq. (|46|l we have to calculate the 
three pre-exponential factors E+(t) — E_(t), [a+(r) — 
Qi_(r)]/d+(r) and o,~ +1 (t) in Eq. ijlfi}. We assume that 
the Hamiltonian H(t) is analytic in r. Let r c be one of 
the branch points of E + (t) — E-(t). Close to r c we get: 

E+(t) - E_(t) = c(r - rc) 1 ' 2 (84) 

with c a constant, depending on t c . Eq. (|51|) implies 

z-z c = ^ c (r-r c ) 3 / 2 , (85) 

where z c = z(t c ). In the adiabatic limit S — > the 
main contribution to the integral (2. line of Eq. (|46l) ) 
comes from the singular points r c and z c , respectively. 
Consequently we have to calculate the pre-exponential 
factors (2. line of Eq. 1)461) ) close to the singularities, 
only. Let us start with [oq_(Y) — o:_(t)]/&+(t). Using 
Eq. H50[) it follows with a±(z) = a±(r(z)) close to z c 

[a+{z) - a-{z)]/a' + {z) S 6(z - z c ), (86) 
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where ' denotes derivative with respect to z. Note that 
H\\,Hi2 and H 2 2 do not enter in Eq. (|86[l . The cal- 
culation of a~ +1 (r) is more evolved. As a first step we 
eliminate 6 m _ 1 (r), b m _ 1 (r) and 6^(r) from Eqs. (J2HJ, 
(12411 which yields a recursion relation for a~ (r) : 



E+(t)-E-(t) 



K_(r) 



K_(r)K + (r)a m _ 1 (r)] 



(87) 



Next we calculate the various terms close to t c . From 
Eqs. (jUnj and lO we get: 



Mr) 
«_(t) 



and 



k_(t)k + (t) = ^( T - r c) 2 - 



(89) 



Expressing the r-derivatives of a m and ct m _i by deriva- 
tives with respect to z: 



a m ( T ) - c{t -T c ) 1/2 a m {z) 



(90) 



and 



a m (T)^c 2 {T-T c )a m '(z) + -(T-T c ) 1/2 a m (z), (91) 

where dz/dr = E + (t) — E-{t) and 1|84[) was used, we get 
from Eq. (E3 with Eq. JBEJl, JEJ3>: 



a m 0) - (-*) 



a m-l( z ) a m-l( 2 ) 



36(z - Z c ) 2 (z - Z c ) 

with initial condition (cf. Eq. H28|) h 
flg (z) = 1. 



-»m-l(*) 



(92) 
(93) 



The recursion relation is identical to Eq. I|30[) in Ref. , 
except the different sign in front of the square bracket. 
The sign change is irrelevant. The exact solution of 
Eq. {3H can be taken from Ref. p^ : 



with 



a m (z) = B m (z - z c y 



(m- l)\(m- I)! 



(94) 



m!(-|)!(-|)! ■ 



(95) 



Then we get from Eq. ggj) with {M jl - (fHH ft and {SSJ 



<*r «n + i(T) a " (T i "~ (T) [^ + (r) - B-(r)]/(r) 
a + (r) 

_/ a + (z) - a_(z) 
<fo a ™+l( z ) TTTTa 

c 

= -6(n + l)B n+1 J2jdz {z _^ k))n+1 I(z) 



= -6(n + l)27«(i) B n+1 J2Hzc{k)), 



(96) 



where the sum over k is restricted to all singular points 
z c (k) above the contour C = {z = z{t)\ — oo < r < oo}. 
Since 



6 • 27r(n + l)B n+l /n\ — > i n+ \ n ^ oo, (97) 



the prefactor in front of m Eq. (|96l) equals (— 1)" +1 . 

fc 

Substituting into Eq. l(4^|) 5 n cancels and one obtains 
the result Eq. (52)l. 
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